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ABSTRACT 



This study addresses a general class of decision situations whose solutions are 
directly applicable to military inventory acquisitions and or disposals. Although 
optimal solutions are well known when subsequent costs are linear to the amount 
of surplus or shortage, the perhaps more realistic case of non-linear costs has not 
been extensively studied. The results of this study suggest the optimal solutions, 
i.e., acquisition quantity or supply, for both conditions of risk and uncertainty about 
demand when the associated cost function is non-linear and quadratic. For 
conditions of risk, optimal solutions are found which will yield minimum expected 
costs for the two-piece cost function where surplus and shortage costs are 
quadratic. This is done for both discrete and continuous demand variables. When 
future need for the item is unknown and only the maximum value can be estimated, 
optimal solutions are obtained for goals of minimaxing cost, minimaxing regret, 
and the Laplace criteria using a uniform probability distribution. It is shown that 
these different approaches to determining acquisition quantities under conditions 
of uncertainty lead, for this general class of decision problems, to the same optimal 
result. Hopefully, this information will aid in the decision process while making 
affordability assessments of new acquisitions. 
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EXECUTIVE SUMMARY 



This study addresses a general class of decision situations whose solutions are 
directly applicable to inventory acquisitions and or disposals. Although optimal 
solutions are well known when subsequent costs are linear to the amount of surplus 
or shortage, the perhaps more realistic case of non-linear costs has not been 
extensively studied. 

Even the most casual observer of the defense marketplace would agree that 
long-term growth in the costs of military hardware has been substantial. Agencies, 
such as the Center for Naval Analysis, are often involved with making affordability 
assessments of new acquisitions under the Department of Defense. In determining 
the affordability of a new acquisition, the entire life-cycle costs need to be 
evaluated. When acquiring a major piece of military hardware, costs associated 
with items such as spares or replacement parts can consume a significant portion 
of the total costs. The timing of when resources are expended plays a crucial role 
in determining the affordability of a new program. Spares or replacement parts not 
purchased at acquisition time but needed (and purchased) later, tend to cost much, 
much more. Furthermore, when the particular piece of hardware is at the forefront 
of new technology the reliability of it or its components may be unknown. 
Estimates of the probability distribution of demand for the number of replacement 
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parts over the life of the hardware will have to be considered along with the 
associated cost functions. 

The results of this study suggest the optimal solutions, i.e., acquisition 
quantity or supply, for both conditions of risk and uncertainty about demand when 
the associated cost function is non-linear and quadratic. This is done for two-piece 
cost functions where each part (surplus and shortage) recognizes both linear and 
quadratic terms. The optimal solutions found for conditions of risk will yield the 
minimum expected costs. For the continuous case, it is derived by differentiating 
the expected cost equation with respect to the decision variable, i.e., the acquisition 
quantity, and setting the result equal to zero. For the discrete case, it is found that 
sufficiency conditions for the local minimum of a function (finite difference 
equations) are elemental in deriving the optimal acquisition quantity. When future 
need for the item is unknown and only the maximum value can be estimated, 
optimal solutions are obtained for goals of minimaxing cost, minimaxing regret, 
and the Laplace criteria using a uniform probability distribution. 

This analysis of a quadratic cost equation has not been studied extensively in 
the past. Hopefully, this information will aid in the decision process while making 
affordability assessments of new acquisitions. 
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I. INTRODUCTION 



Stochastic single-period inventory models, more commonly known as 
"Newspaper Boy" problems, first made an appearance in 1950. where they were 
described in Morse and Kimball’s early Operations Research text [Ref. 1], Based 
on certain assumptions about demand, the problem involves making inventory 
acquisition decisions with the objective to minimize expected cost. These models 
are directly applicable to determining the optimal quantity of spares to acquire with 
a new system [Ref. 2], 

When acquiring a major piece of military hardware, costs associated with 
spares or replacement parts are a significant portion of the total costs. Replacement 
items not purchased at acquisition time but needed later, tend to cost much, much 
more. Furthermore, when the particular piece of hardware is at the forefront of 
new technology the reliability of it or its components may be unknown. Estimates 
of the probability distribution of demand for the number of replacement parts over 
the life of the hardware will have to be considered along with the associated cost 
functions. Estimating the demand uncertainty was researched in a study by Pranom 
Srinopakoon [Ref. 3]. His study suggest some procedures which can be used to 
estimate the demand distribution even when data on unsatisfied demands are not 
available. In another study, a decision procedure for the transition from uncertainty 
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to risk was proposed by Kadir Sagdic [Ref.4] who explores the use of order- 
statistic-based quantile estimators as a decision procedure while data accumulates 
in the early decision periods. 

In existing applications, cost functions have been assumed to be linearly 
proportional to quantity. Although optimal solutions are well known when 
subsequent costs are linear to the amount of surplus or shortage the perhaps more 
realistic case of non-linear costs has not been extensively studied. Little research 
has been done when associated costs are non-linear. However, even the most 
casual observer of the defense marketplace would agree that long-term growth in 
the cost of military hardware has been, and maybe expected to be, substantial. The 
purpose of the study reported in this thesis is to seek and examine solutions where 
costs behave in a non-linear manner, i.e., quadratically, as functions of shortage or 
surplus. 

Chapter II briefly reviews the standard model where costs are linearly 
proportional to surplus or shortage. Chapter III is devoted to the purpose of this 
thesis, i.e., to the quadratic extension of the newsboy problem. Hence, optimal 
acquisition decisions will be sought for both conditions of risk and uncertainty 
about demand. Three principles of choice will be applied to decisions under 
uncertainty: the Laplace principle, minimax cost solutions, and minimax regret 
solutions. Conclusions and recommendations are given in the final chapter. 
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II. SINGLE-PERIOD INVENTORY MODEL (LINEAR COST 
FUNCTIONS) 

In this section we shall briefly describe the well-known newsboy problem, 
illustrating its optimal solutions under both risk and uncertainty with respect to 
demand. This generalized problem appears frequently in a variety of scenarios. 
The model may represent the inventory of an item that (1) becomes obsolete 
quickly, such as the daily newspaper; (2) spoils quickly, such as vegetables or 
loaves of bread; (3) is stocked oniy once, such as spare parts for a single 
production run of a new piece of military hardware; or (4) has a future that is 
uncertain beyond a single period, such as designer clothing fads [Ref. 5]. 

In general the inventory model has been structured such that orders are placed 
once for purchase of inventory to cover a single period. Stock shortages may not 
be refilled and stock surplus may not be transferred for use in the next period. 
Costs incurred by shortage or by surplus are proportional to the difference between 
the quantity on hand and the subsequent demand during the period. If exact 
demand were known it would be a trivial matter to minimize cost: one would 
simply acquire that exact amount. Unfortunately, since demand is typically 
unknown at the time of acquisition, and could be characterized as a random 
variable, then minimizing expected cost becomes a more complicated task. 



3 



However, for an appropriate decision rule, an optimal inventory policy can be 
obtained even if limited information about the probability distribution of demand 
is available. It is clear that ordering some amount more than the possible minimum 
demand would be desirable, but certainly the order should be no more than the 
maximum demand. 

The single-period cost equation may be structured from the following 
components. 

Q = Quantity (supply) of inventory initially to be on hand 

for the period. This is the decision variable for the problem. 

D = Demand for the inventory during the period. Demand is 

unknown at the time supply Q is selected, and may usefully be 
considered as a discrete or continuous random variable. 

C s = Cost of surplus for each unit of inventory unsold, i.e., when 
supply Q exceeds demand D. 

C 0 = Cost of outage for each unit of unsatisfied demand, i.e., when 
demand D exceeds supply Q. 

p(d) = Probability distribution of demand, if demand is discrete. 

f(d) = Probability density function of demand, if demand is continuous. 



The cost function, depending on whether there is a surplus or a shortage is 



Cost=C[Q) = 



' C S (Q-D) , 
C Q ( D-Q ) , 



OS D <>Q, 
0<D. 



( 1 ) 



4 



For simplicity in familiarization with the basic model we’ll assume that both Q and 
D are continuous variables. The discrete version will be explained in the next 
chapter, dealing with the non-linear cost function. The linear cost function appears 
graphically as in Figure 1. 




Figure 1. The Newsboy Cost Function 



Given a certain density function f(d) for the continuous random variable of 
demand, and seeking a value for Q\ the decision variable that the decision maker 
is trying to find over the range between zero and D max , the expected cost is 
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E[C(Q)] = j (Cost) f(D) dD. 

0 

This, upon including the cost function (1) becomes, 

oo oo 

E[C(Q)] = jc a (Q-D) f (D) dD + fc o (D-Q) f(D) dD. (2) 

o o 

We seek the value Q\ which minimizes this expected cost function (2). The 
optimal value of Q may be found by differentiating the expected cost with respect 
to Q, equating the derivative to zero, and solving for the optimum value Q*. Since 
Q appears in the limits we may use Leibnitz rule [Ref. 8] for integration, or 

J = Jc g f{D)dD + J-C Q f(D) dD = 0 . 

o 

Applying the cumulative distribution function 

X 

Fix) = j fit) dt 

— OO 

to obtain 

c a FiQ*) - C 0 [l - FiO *) ] = 0, 



we arrive at the well-known result [Ref. 5] 
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F(Q') 



( 3 ) 



C. + C„ 



To verify that the expected value of cost is minimized at Q‘ we should 
examine the second derivative 



d 2 g[C(g) 1 

dO 2 



o-o’ 



(C a + C Q ) f(0*) • 



Since the quantity (C s +C 0 )f(Q*) is positive the value indeed is minimized. 

To summarize, the above result says that in order to minimize the expected 
value of cost, the decision maker should order the quantity Q* which is equal to the 
[C^Q+CJJth quantile of the demand distribution, as shown in Figure 2. [Ref. 6] 

Furthermore, as a decision tool, it is not necessary to know the exact 
numerical values of C s and C 0 , since the ratio of the two may be sufficient. For 
example, suppose that as a best cost estimate a unit of surplus will cost twice as 
much as a unit of shortage, i.e., C s = 2C 0 , which equates to determining the value 
for the .33 quantile estimate. 



I 
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When the distribution of demand is unknown the decision becomes one under 
uncertainty. Three well-known approaches to decision making under uncertainty 
are the Laplace solution under the assumption of a uniform distribution of demand, 
minimax cost solutions and minimax regret solutions. For the newsboy problem, 
these all lead to the same rule to find the optimal solution. We shall discuss only 
the minimax cost approach. [Ref. 8] 

This approach is to choose Q so that the worst possible cost will be 
minimized, that is, to minimize maximum cost. As is evident from the solid line 
in Figure 3 the maximum cost occurs at D = D max . By increasing Q to say Q’ we 
can reduce the cost at D = D max , as shown by the dashed line in Figure 3. Note, 
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however, that the maximum cost now occurs when demand is at D=() and so the 



process continues. 




Figure 3. Illustration Of Extreme Point Maximums Of The Newsboy 
Cost Function With Bounded Demand. 

Clearly, we would want to choose Q so that the cost at D = D max is equal to the 
cost at D = 0. Equating these, we obtain 

- Q*) 
or 




9 



as the order quantity Q* which minimaxes cost. 

In this chapter we have familiarized ourselves with the minimum expected 
cost solutions under both risk and uncertainty conditions of the basic newsboy 
problem. We will now proceed to the more extensive case where associated costs 
are non-linearly proportional to the amount of surplus or shortage on hand. This 
will be done for both continuous and discrete quantities of supply and demand. 
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III. SINGLE-PERIOD INVENTORY MODEL (QUADRATIC COST 
FUNCTIONS) 

This chapter is devoted to the quadratic extension of the standard single- 
period model, i.e., the newsboy problem. Costs are no longer assumed to be 
linearly proportional to the difference between the quantity on hand at the 
beginning of the period and the subsequent amount demanded throughout the 
period. Instead, costs are assumed to increase quadratically in relation to a surplus 
or shortage. With the exception of this new cost relation the other characteristics 
of the standard model remain unchanged; i.e., orders are placed once for the 
purchase of inventory to cover a single period, stock shortages may not be refilled 
and stock surpluses may not be transferred for use in the next period. The cost 
function may be written as two second degree polynomials in the form of 




0 <LD±Q 
Q<D, 



( 5 ) 



for the continuous case, and 




for the discrete case. 
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The components for the quadratic cost equation are as follows: 

Q = Quantity (supply) of inventory initially on hand for the period. 
This is the decision variable for the problem. 

D = Demand for the inventory during the period. Demand is 

unknown at the time supply Q is selected, and may be usefully 
considered as a discrete or continuous random variable. 

C<, = Cost of surplus for each squared unit of inventory unsold, i.e., 
when supply Q exceeds demand D. 

C s2 = Cost of surplus for each single unit of inventory unsold. 

C ol = Cost of outage for each squared unit of unsatisfied demand, i.e., 
when demand D exceeds supply Q. 

C o2 = Cost of outage for each single unit of unsatisfied demand. 

The cost coefficients may assume both positive and negative values. 
However our work will address cases where the four cost coefficients are non- 
negative. 1 Under these conditions of a quadratic cost equation the penalty of 
having a surplus or shortage is much greater than it is when the cost equation is 
linear. Total cost increases rapidly when the quantity ordered does not coincide 
with the subsequent demand. Notice that the cost function is a strictly convex 
function and may appear as in Figure 4. 



'Cases where C s2 =C o2 =0 have been investigated [Ref. 10, 11]. 
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Figure 4. Illustration Of The Quadratic Cost Function 
A. SOLUTIONS UNDER RISK 

Our ability to find good solutions for the optimal order quantity depends upon 
the amount of information we have about the future; in particular, about the 
magnitude of demand D. In this section we will assume that the probability 
distribution of demand is known. Our goal is to find the value of the decision 
variable which minimizes expected cost. 

1. Minimizing Expected Cost, Continuous Case 

Given that demand D is a continuous random variable with a known 
density function f(d), and that Q (the decision variable for the optimal order 
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quantity) is also a continuous variable; the expected cost equation is 



E{C(Q)] = | (Cost) f(D) dD. 

D= 0 



Which, upon including the cost function (5) becomes 



o 

E[C(Q) ] = f C Si ( Q-D ) 2 f (D) dD 

D= 0 
0 

+ f C Si (Q-D) f (D) dD 

D= 0 

CD 

+ f C 0 i (D-Q ) 2 f(D) dD 

D=Q 

m 

* f C 0 j (D-Q) f (D) dD. 

D-0 



( 7 ) 



As with the standard single-period model, the optimal value of Q may be found by 
differentiating the expected cost equation with respect to Q and then equating the 
derivative to zero. Since Q appears in the limits of integration, we may use 
Liebnitz’ rule for differentiation of integrals [Ref 7]. The result is 
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dE[C(Q) ] 

dQ 



o 

= f 2 C Si (Q-D) f(D) dB + C Si (Q-0) 2 f ( 0) 

D=Q 

0 

+ j f (D) dD + C S2 (Q-0) f(0) 

D=0 

+ f (-2) C Qi ( D-0 ) f(D)dD + C Qi (Q~0) 2 f(D) 

D=0 

co 

+ [ -C 0i f(D) dD + ( 0-0 ) f(D) = 0 . 

0=0 



Which further simplifies to 



E[C(Q) ] = 

dQ 



2 C Si 0F(0) ~2C Si 

+ C„ F(Q) 
*2 



Q 

f D f(D) dD 
o=o 



+ 2 C 0i 0 [1 ~F(0)] ~ 2 C Qi J D f (D) dD 

- [1 -F(0)] =0. D ° 



By adding and subtracting the quantity 

o 

2 C 0i f D f(D) dD, 
o=o 



the final result is 
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( 8 ) 



F(Q) [2Q{C Sx -C 0 ) +C 3 ^C 0 ] + 2C^{Q-E[D] ) 

Q 

- C 02 - 2 (C St - C 0i ) f D f(D) dD - 0 . 

D= 0 



The expected cost is minimized for the value of Q = Q* which satisfies this 
condition. 

Depending upon the distribution, solving (8) for Q* may be a tedious 
task accomplished best by a series of iterations. On the other hand, the distribution 
may allow us to be able to solve for Q* by simple manipulation. As an example, 
suppose that an elemental replacement part is available for a new piece of military 
hardware where the mark-up for this part is relatively high and the penalty for 
being out of it is unforgiven by the consumer. Let’s suppose that each unit of 
shortage (in thousands of dollars) costs eight for a single unit with an additional 
cost of two for missing the subsequent demand, i.e., C 0l = 2.00 and C o2 = 8.00. 
Each unit of surplus (in thousands of dollars) costs one and its additional cost is 
ten percent for missing the subsequent demand, i.e., C sl = 0.10 and C s2 = 1.00. 
Thus the cost equation is 



Cost - C(Q) 



0.10 ( Q-D f + 1.00 ( Q-D ), D=0,l,...Q 

2.00 (. D-Q f + 8.00 (D-Q), £=<?+!,.... 



We’ll assume, that the item has a demand that is exponentially distributed with a 
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mean of 200, i.e., E[D] = 200. By substituting the values for the coefficients and 
applying the following two characteristics of the exponential distribution [Ref. 1 1], 

y_ 

F(Q) = 1 - « 200 
and 




e 200 dy = 200 - Qe 200 



__Q_ 

- 200 e 200 



to Equation (8) the result becomes 



_Q_ 

[: l-e 200 ][2Q(-L9)+9] + 4(a-200) 

__Q_ Q 

- 8 - 2(-1.9 )[-Qe 200 -200e 200 +200] = 0. 



When simplified, this expression reduces to the following equality: 



__Q_ 

0.20 Q - 769 e 200 = 39 . 



As we can see from Table 1, the value for Q which satisfies this 
equality lies between 504 and 505 units. Obviously, under the above constraints, 
it is less costly to have a surplus than it is to have a shortage. However, it is 
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interesting to note that here the optimal order quantity is more than two and a half 
times the expected demand quantity. 



Table I. Optimal Order Quantity For An 
Exponential Demand Distribution 



Q 


. 2Q - 7 


69e -Q/2°° 


200 


-242 


1.89 


400 


-24 


.07 


450 


8 . 


94 


500 


36 . 


87 


500 


38 . 


92 


505 


39 . 


43 


510 


41. 


95 


550 


60 . 


83 



In comparison, if there are no costs associated with the quadratic terms, i.e., the 
model is a simple "Newspaper Boy" problem which has a linear cost function, then 
the optimal order quantity would be approximately 440 units. This value is also 
much greater than the value for the expected demand; however, it is somewhat less 
than the optimal value for the quadratic cost function. The quadratic solution 
requires a larger supply quantity because the costs associated with a shortage are 
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more substantial than the costs associated with a surplus. 



2. Minimizing Expected Cost, Discrete Case 

In most situations inventory exchanges are measured in discrete units 
where quantities are too small to permit approximation by continuous variables. 
In this case, both Q and D are treated as discrete variables, where p(D) is the 
probability distribution of demand. The expected cost for ordering Q items 
E[C(Q)J, is 

ao 

£[C«?)] - £ (Cost) p(d). 

D= 0 



This, upon including the discrete cost Equation (6), becomes 



£tC(«] - £ (C.(Q-Df p(D) * CAQ-D) p(D)) * 

o-o (9) 

OB 

E [C. (»-<?)’ ptf>) ♦ C (D-Q)p(D)]. 

D-Q+l 



In terms of the expected cost function, sufficient conditions for a local minimum 
of a function for Q = Q* and defined for integer values of its argument are 

£[C«T + 1)] - £[C«?*)] ^ 0 ( 10 ) 

and 
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(11) 



The forward difference. Equation ( 10), and the backward difference. Equation (11), 
will be used to construct the decision rule to minimize expected cost. 

Beginning with the forward difference, the expected cost, using Equation 
(9), at Q+l is 



<?+ 1 

E[C(Q+ 1)] = £ C s (Q+l-D) 2 p(D) + 

D=0 1 

Q* 1 

£ C AQ+l-D) p(D) + 

/)= 0 
00 

E C (D-<?-l) 2 /-(£>) + 

D=Q+2 

oa 

E C (C-0-1) /.(D) . 

D-Q+2 



When expanded this becomes. 



<3+1 <3+1 <?+i 

£[C«?+1)) - E C, (.Q-D) 2 p(P) * E 2C (Q-D)p(D) +£C p(D) 

D=0 ‘ D=Q D = 0 

Q 

+ c «?+l-<?-l) P(D) + £ C (Q+l-D) p{D) + 

jD=0 



E C <D-Q?p(P) - E 

D=<3+2 D=Q+2 

CO 

+ E C (D-® p(D) - 
£>=<?+ 2 



2C 0 (D-Q)p(D) 

CO 

E C„p( 2 >) 

Z)*0+2 



+ E 

D=Q+2 
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If the expected cost at Q is truly a minimum, then the difference between the 
expected cost at Q+l and the expected cost at Q must be non-negative. Therefore, 
upon applying Equation (10), the above result (after some effort) becomes 



Q 

2 tE CAQ-D) p{D) - £ CAD- 

D=0 



D=<?+ 1 



AC. * CJE p(D) + (C - C ) £ KD) *0 

D=Q+1 



D= 0 



By using the well known identity 



«<?) - Ep(°) 

D = 0 



1 - £ />(£»), 

£>=<?+! 



the equation simplifies to 

2[E c iQ-mm- E c (D-® P (i>)] 

D =0 D =<?+1 ( 12 ) 

i(- C,, - + C ci - C c )P(Q) * (C^-Cy. 

Equation (12) satisfies the first half of the sufficient conditions for the construction 
of the decision rule. 

Next, we’ll apply the backward difference. Using Equation (9) again, 
the expected cost at Q-l is 
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CH 



E[C(Q- 1)] = £ C, (Q-l-D) 2 p(D ) + 

D=0 ' 



CM 



]T C «?-l-D) p(D) + 

D=Q 



£ C (D-Q-1) J />(D) + 

D-Q 



£ C (D-O+l) p(D) . 

D-Q 



When expanded and simplified this becomes, 



Q Q Q - 1 

£[C«?-1)] - £ C,<Q-D? PV>) -2 £ CJQ-DMD) + £ C,jm 



D= 0 



D=0 



D=0 



Q - 1 



+ E - E p( fl ) + 

D=0 D=0 



£ C (D-Q) z p(D) +2£ C (D-Q)p(D) + £ C P (D) 

Z>=<? /)=(? 2)=<? 

OO OD 

+ E CJD-Q)p(D ) * £ C ->(£>). 

D=<? D-Q 



As before, if the expected cost at Q is truly a minimum, then the difference 
between the expected cost at Q-l and the expected cost at Q must be positive. 
Applying Equation (11), the above result becomes 
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Q-l 



S, Sj- 






D = 0 



°l °2- 



D=Q 



21E CAQ-D) p(D) - Y C ( 

D=0 D=Q+1 



By using another well known identity 



-1) = f>(D) = 1 - f>(D), 

D=C D=Q 



the equation further simplifies to 



(c,, - c Si - c 0| - cyp«?-i) - (c^cy 

Q 09 

> 2 [£ CAQ-D) p(D) - Y CAD-Q) p(D )]. 

D=0 D=Q+l 



(13) 



Equation (13) satisfies the second half of the sufficient conditions for the 
construction of the decision rule. Finally, by combining Equations (12) and (13), 
the first order difference conditions, we arrive at the decision rule 

(C,, - c, 2 - C 0i - cyw-i) + (C„_ + cy > 

Q » 

2tE C AQ-D)p(m- Y c o (d-q) p(d)] k < |4 > 

D=0 1 0=<M 1 

(- c - C + C - C )P(Q ) + (C - C ) . 
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The optimal value of Q = Q‘ which satisfies Equation (14) will yield the minimum 
expected cost solution to Equation (9). 

Finding the value of Q which satisfies this double inequality may be a 
rather intractable task. In some cases, rather than use (14) it may be easier to 
handle the discrete case by direct enumeration of expected costs for each possible 
value of order quantity Q. For example, suppose that C sl = 2, C s2 = 4, C ol = 3, C o2 
= 6, and that the probability distribution of demand and cumulative probability 
function of demand are, respectively. 



0.1 


D= 0, 






0.1 


D= 0, 


0.2 


D= 1, 






0.3 


D= 1, 


0.4 


D =2, 


and 


P(D)= - 


0.7 


D= 2, 


0.2 


D= 3, 






0.9 


D= 3, 


0.1 


£>= 4, 






1.0 


£5=4. 



Values of the cost function associated with the various values of supply Q, and 
demand D, are as shown in Table II. For example, if supply Q, at the beginning 
of the period, is equal to 3 and demand D, by the end of the period, was only 1, 
the cost would be 

CostiQJ)) = C(3,l) = 2 (3 - 1) 2 + 4(3 - 1) = 16. 
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Table II. Cost Matrix For A Discrete Quadratic Cost Function. 



DEMAND D 


S 

U 

P 

P 

L 

Y 

Q 




1 

0 


| 

1 


2 


3 


4 


0 


6 


9 


24 


45 


72 


1 


6 


0 


9 


24 


45 


2 


16 


6 


0 


9 


24 


3 


30 


16 


6 


9 


9 


4 


48 


30 


16 


6 


0 



Notice that the cost matrix is an upper and lower semi-symmetric matrix with zeros 
on the diagonal. The expected cost E[C(Q)] may be computed for each of the five 



candidate values of Q as 

£[C(0)]= 0(.l) + 
£[C(1)]= 6(.l) + 
£[C( 2)]= 16(A) + 
£[C( 3)]= 30(.l) + 



9(2) + 24(A) + 45(2) + 
0(.2) + 9(.4) + 24(.2) + 
6(.2) + 0(.4) + 9(2) + 
16(.2) + 6(.4) + 0(.2) + 
+ 30(.2) + 16(.4) + 6(2) + 



72(.l) = 27.60 
45(.l) = 13.50 
24(.l) = 7.00 
9(.l) = 9.50 
Q(.l) = 18.40. 

optimal supply value is 



E[C(4)]= 48(.l) 



We conclude that, by direct enumeration, the 
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Q* = 2. Notice that, even though the cost for having a shortage is significantly 
greater than it is for having a surplus, in this case the optimal supply value is the 
same as the expected demand value. If Equation (X) were used to solve for the 
optimal value Q then the task would be more cumbersome. However, it should 
also lead to the same solution for the optimal supply value Q*. For example, by 
substituting the values for the cost coefficients. Equation (14) becomes 

Q 00 

-11 P(Q- 1)+9 > 4^ (Q-D)p(D ) -6 £ (D-Q)p(D) z -9 P(Q)~ 3. 

D =0 D =<?+1 

As we can see from Table III, the optimal supply value Q that satisfies Equation 
(14) is Q* = 2, which is the same as before. 
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Table III. Optimal Supply Value Q* For A Discrete Cost Equation. 



Q 


-IIP (Q+l) +9 


4 I(Q-D)p(D)-6S(D-Q)p(d) 


-9P(Q)-3 


0 


9.0 


4 ( 0 . 0 ) -6 ( 2 . 0 ) = -12 


-3 . 9 


1 


7.9 


4 ( 0 . 1 ) - 6 ( 1 . 1 ) = - 6.2 


- 5.7 


2 


5.7 


4 ( 0 . 4) -6 ( 0 . 4 ) = - 0.8 


- 9.3 


3 


1.3 


4 ( 1 . 1 ) - 6 ( 0 . 1 ) = 3.8 


- 11.1 


4 


- 0.9 


4 ( 2 . 0 ) - 6 ( 0 . 0 ) = 8.0 


- 12 . 0 


5 


- 2.0 


4 ( 3 . 0 )— 6 ( 0 . 0 ) = 12.0 


12 . 0 
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B. SOLUTIONS UNDER UNCERTAINTY 



Deciding on Q when the probabilities of future demand are unknown and not 
estimated falls into the class called decisions under uncertainty . Here, although the 
probabilities are unknown, it may be possible to have some information about the 
distribution without actually knowing the distribution. For example, we may be 
able to estimate one or more descriptive values relating to the demand distribution, 
such as mean, mode, or median. We may have an idea about the shape of the 
distribution as being symmetric or perhaps skewed. We may also be able to 
estimate the bounds for the range of demand values; e.g., demand could range from 
zero to say D max . All of these provide information which may be useful in 
determining the optimal supply value. 

We shall examine the last case, where the only estimate that can be made 
about the nature of the unknown demand is that its range is from zero to some 
upper bound D max : . 



^he general case of D min < D <D max is readily handled by noting that we would 
always stock D mm , and the decision problem concerns the amount of supply in 
excess of D min . 
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Given this as our only information, three types of rationale (principles of choice) 
have been proposed for decisions under uncertainty [Ref. 8], They are 

• The Laplace Solution, 

• Minimax Cost Solutions, and 

• Minimax Regret Solutions. 



Although different in their objectives, we will show that all three lead to the same 
rule for computing Q*. 

1. The Laplace Solutions 

The Laplace approach to the case of uncertainty is to assume that 
demand is uniformly distributed over the demand interval and then use an expected 
value solution. In our case, the interval is 0 < D < D max . We seek to examine both 
cases of continuous and discrete demand. 

a. Continuous case 

If demand is a continuous random variable, we would assume that 
the distribution is rectangular with density 




0 <. B <. D_„. 

TTtflT 
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Thus, the distribution function is 



F(D) = 



D 



max 



Applying this and the identity. 



Q Q 

I D fi4)dD = J D 
£>= 0 2)=0 








to the expected cost solution (8) 

F{Q){2Q(C S -C 0 )+C S +C 0 ) + 2C 0i (Q-E[D]) 

Q 

- C 0i - 2 (C s - C 0 ) f DflP)dD = 0 , 

D = 0 



( 8 ) 



we have 



D. 



+ 2 C g ( 



max 



- C - 2(C - CJ- 

°i K s i °i 2 



max 
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When expanded this becomes. 



(<e*-cj\ 



\ max / 



((C^CA _ \ 



+2C 



ol 



\ max 



" ( C ol D m*x + C o2 ) = 0 5 



which, as we notice, is a quadratic equation in Q. This leads us to 3 



(C^C. 

— — — +2 C 



Q = 



ol 






max 



\ 


/ 


+ 




h 





c s2 +C o2 



\2 



+2C 



ol 



max 



+ 4 



( C -C \ 

^ si ^ ol 



max 



(C 7 D + C 

V ol max o2' 



2 



( C -C N 



\ max 7 



(15) 



As a check, we need to ensure that the rules of mathematics are not violated, i.e., 
the value under the radical must be non-negative and the divisor of a fraction must 
not be zero. By expanding the value under the radical we get 



(C,2^2) 2 . “C OI (C j2+ CJ + 4C ^ 



max 7 



D. 



4C ,C n ,D 4 C ,C , 

5/ 01 max sl o2 



4C 2 ,D 4 C ,C , 

0 / max ol o2 



D. 



max 



D. 



max 



D. 



max 



D 



max 



3 The general solution for quadratic equations provides two solutions, but we are 
only concerned with the non-negative values of Q. 
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which equals. 



^ s2 + 2 ^s2^o2 



C 



o2 



4 ^of'sl 



^ofo2 



i 1 

max 



j 2 

max 



^Ll 



4Q + 4C sl C 01 + 



" 4C 0l 



^C o 1 C o2 



max 



Notice that the last two terms cancel are cancelled by other terms and the 
remaining terms are all positive: therefore, the sign rule under the radical has not 
been violated. For the other rule, if C sl equals C 0l then the divisor of the fraction 
is zero. However, if this were the case we wouldn’t have a quadratic equation in 
Q. The solution for the optimal value of Q = Q‘ would simply be 



Q ' = 



C j£) + C 5 

or mb o2 



' (P*+CJ 
D „ 



+ 2 C 



Ol 



In comparison, this provides a simple adjustment to the Laplace solution of a linear 
cost function which is 



Q * 



( ^02 
(C s2 + C o2 




5 



and where the coefficients of the quadratic terms are equal. 
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In conclusion, for the Laplace criterion. Equation (15) above provides the general 
means of computing the optimal supply quantity Q = Q* which will yield the 
minimum expected cost for a quadratic cost function when a decision under 
uncertainty is made. 

b. Discrete case 

If demand is a discrete random variable, the Laplace approach 
would assume the uniform distribution for demand 



m 




'= 0 , 1 ,..., 



max 



Here, the distribution function is 



PiD) 




+1 

+ 



and for this case the optimal expected value rule (14) 

(C,, - - c Ci - cjPiQ- 1) - (C 0i - cy > 

Q OO 

2[E c , «?-■£» p(P) - E cjp-q) P m 2 

D=0 1 1 

(- ^ - c h + c„ x - c 02 )P(Q) * (C fl2 - cy , 



(14) 
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leads to 



(C - C - C - C )- 

v s i s i °i V 



muT 



+ 1 



+ (C + C ) > 






<r+ 1 N 

\~^max + ^ } 






'"I \ 

max 



2(I> +D 



<?* + ! 



(-C . - C + C - c )— — — + (C - C ) . 
v si s 2 ®i Vn +1 1 °2 V 



max 



2 

^ / 



(16) 



We have used the identities 



£ D p(D ) = £[D] - £ 0 



£»=<?+! 



0=0 



and 



<?* 



E D pw - 

D=0 



1 <?* «?* + !) 

1} +1 2 

max 



to arrive at Equation (16). Notice that the above expression (16) can be further 
simplified to 
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(C - C - C - CJ— ^ 

s i ** °i °2 n +1 

m ax 



(C fl . - CJ > 






o*+ 1 N 



J) +1 

\ max y 



+ 2 C oI Q--C„ D ^ i 



(-C, - c, + C - CJ-& + (C - CJ , 
v Ji s 2 °t Yn +i v °2 V * 

max 



and ultimately to 



^sjQ + C s2 + c sI > 

CJD^- Q*) 2 + CJP^- <?*) - C sJ Q * 2 - C s2 0* * 

2 C^OU- <T) - C oJ - C o2 . 



( 17 ) 



This is the decision rule which leads to the optimal value Q* which must 
simultaneously satisfy two quadratic expressions in Q, i.e., 

(Csl ~C 0 i)Q* 2+ @C sl + C s2 + C o 2 )Q+(C s1 +C s 2 - ~C 0 2®m*x ) * 0 

and 

(C S 1 -C 0 j)Q * 2+ (^ s2 + 2C oi D nuix +2C oi + C o2 )Q-[C ol (D taax + l) 2+ ^ 0 2(^max + l)5 > ^* 



This value for which Q=Q* will yield the minimum expected cost for the quadratic 
cost function when a Laplace decision under uncertainty is made and the only 
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information we have about demand is that it is less than or equal to, say, D max . 



2. Minimax Cost Solution 

Another approach to decisions under uncertainty about future demand 
is to choose Q so that the worst possible future cost will be as small a possible, 
that is, to minimize maximum cost. 

a. Continuous case 

For the case of continuous demand, the optimal value of Q in 
terms of minimaxing cost is easily found. As shown graphically in Figure 5 we 
can see that for the quadratic cost function at the indicated value of Q, the 
maximum cost will occur when demand is equal to D max . By increasing Q we can 
reduce the cost at D = D max . If we use a higher value, say Q\ then we will have 
the cost function shown by the dashed line in Figure 6. Note, however, that the 
maximum cost now occurs when demand is at D = 0, and to reduce this maximum 
possible value, we would want to reduce Q from the Q’ value, and so the process 
continues. It should be clear that we will minimize the maximum cost when Q is 
chosen so that the cost at D = 0 is equal to the cost at D = D max . This is the 
approach to the minimax cost solution. 
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Figure 5. Maximum Cost In The Quadratic Cost Function With 
Bounded Demand 




Figure 6. Illustration Of Extreme Point Maximums Of The 
Quadratic Cost Function With Bounded Demand. 
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Equating the cost at the two extreme points of demand, we get 



C s jQ* 2 + C s 2 Q* = + 

Notice that this is a quadratic equation in Q and by grouping the terms we have 

(C,i-Coi)Q' 2 * (C* - 2 S_C d + CJQ - C ol D^ - C 0 fi^ - 0. 



If we divide through by D max (not equal to zero) and solve for Q we have 4 



C +C 
<u j2 +<u o2 



D 



+2C 



ol 



max 



\ 


/ 


+ 






A 



C s2 + ^o2 

D 



\2 



+2C 



ol 



max 



+ 4 



/r -C 

^ si ^ ol 



D 



K C^D m *CJ 



max 



( C ~C ^ 
Wi ^ ol 

D 

\ max / 



(18) 

Notice that Equation (18) is identical with the Laplace solution (15) found earlier. 
This is the optimal value of supply (Q=Q*) that yields the minimum expected cost 
for the quadratic cost function when a decision under uncertainty is made. 



4 As before, the general solution for quadratic equations provides for two 
solutions, but we are only concerned with non-negative values of Q. 
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b. Discrete case 



When demand is discrete, the minimax cost solution is not directly 
apparent. As with solutions under risk, in order for costs to be minimaxed at 
Q = Q\ the following inequalities must hold: 

Maximum C(Q *) < maximum C(Q * + 1), (19) 

and 

Maximum C(Q *) < maximum C(Q * - 1). (^0) 

As we saw in the continuous case, the cost function will have its maximum value 
at either C sl (Q) 2 +C s2 (Q), or at C 0l (D max -Q) 2 +C 02 (D max -Q). These reflect the costs at 
the extreme values of D, i.e., D equal to 0 and D equal to D majl . The first condition 
(19) may be written as 

MAX [ C sl (Q *) 2 +C s2 (Q *) , CJD^-Q-f+CJD^-Q') } < 

MAX [ CJQ-+1)*+CJQ'+1) , CJD^-Q'-lf+CJD^-Q--!) ] . 

(21) 



Because 



CJQ-f + CJQ’) < CJQ-+l) 2 + CJQ-+l) 
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and 






,2^-Q’- « < CoAD^-Q ') 1 * CJD^-Q’), 



C s1 (Q’+1) 2 +C s 2 (Q*+1) is the maximum of the right-hand side of (21). (See 
Appendix A for further clarification.) Therefore, 



< ^((r+i^c^Kr-D, < 22 ) 



provides a useful form of the first condition for a minimax cost solution at Q*. 

Proceeding to the second condition (20), it may be written as 

MAX[C SI (Q-)* + CJQ-) . CJD^-Q'f + CJD^-Q-) } < 
»UX[C„(Q‘-V? + CJQ--l) , C t /D, II -(3M) ! *C, 1 (D„-(?M)} 



(23) 



Because 



C, J (<? , -0) 2 «C^«?--0) > C S] (Q “ - l) z + C^iQ * - 1). 



* i\2 



and 






C 0 i(D max -Q*+l) 2 + C o2 (D max -Q*+l) is the maximum of the right-hand side of (23). 
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(See Appendix A for further clarification.) Thus, we have 



C. 






+ C I2 (0*-0) < (24) 



max 






max 



as the useful form of the second condition for a minimax cost solution at Q*. 

Together, from the first condition (21) and the second condition 

(24), we obtain 



2C„Q' + Crf + C tl > 

CJP m - Q’f + CJD^- Q-) - C sl Q' 2 - C s2 Q * > 
2 CJD^- <?•) - C„ - C o2 , 



(25) 



as the optimal decision rule to minimax costs for decisions under uncertainty with 
discrete demand. Notice that the result (25) is the same as that obtained via the 
Laplace approach (17). 

3. Minimax Regret Solutions 

Yet another way to address decisions is through consideration of regret. 
The minimax regret approach to decisions under uncertainty considers the regret 
the decision maker will feel later if the best decision, i.e.. the one which has the 
least cost, were not made. Thus, regret is the difference between the cost incurred 
by choosing the quantity Q and the lowest possible cost that could have been 
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obtained, had the optimal value Q’ been chosen. Therefore, for each value of 
demand. 



Regret (Q) = C(Q ) - minimum {C{Q)\. 

Q 



However, 



minimum [C(Q)] = 0 

Q 



for all values of demand where 

0 * D * 

and thus 

RegreKQ) = C(Q) . 

It follows that for the quadratic cost function, the solutions which minimax regret 
are the same as those which minimax cost, viz., (18) for continuous demand and 
(25) for discrete demand. 

In summary, although different in their objectives, all three (principles of 
choice lead) lead to the same rule for computing the optimal value of supply Q 
which will minimize the expected cost when a decision under uncertainty is made. 
Thus we have obtained a useful rule for selecting the optimal acquisition quantity 
when future demand is not known. 
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IV. DISCUSSION, CONCLUSION AND RECOMMENDATIONS 



Determination of acquisition quantities is of increasing importance for today’s 
Navy. With a declining defense budget, considerable attention is being focused on 
various programs and capabilities. According to a recent report published by the 
Center for Naval Analysis, affordability assessments are being formalized and 
incorporated into the acquisition management process. While the principles and 
methods applicable to those assessments are still in a rudimentary stage, it is clear 
that affordability, which is concerned with top-line budget constraints and the 
magnitude and timing of competing uses of resources, is a macro concept. [Ref. 12] 
The timing of when resources are expended plays a crucial role in 
determining the affordability of a new program. As stated in the introduction, 
when acquiring a major piece of hardware (such as a tactical aircraft) one usually 
can buy spare parts at the same time. Replacement items not purchased at that 
time but needed later tend to cost much, much more, due to such factors as 
retooling, lost technology, and down time while the item is being acquired. 
Resources allocated to the acquisition of spare parts for a major piece of hardware 
can consume a significant portion of the item’s total budget; therefore, making it 
unaffordable. 
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In a book written by Norman Augustine, Augustine's Laws , he finds a 
colorful way of characterizing the growth that has occurred in the cost of military 
aircraft: 

In the year 2054, the entire defense budget will purchase just one aircraft. 
This aircraft will have to be shared by the Air Force and the Navy 3-1/2 days 
each per week except for leap year, when it will be made available to the 
Marines for the extra day. [Ref. 13:p.l43] 

The decision of the optimal order quantity of spares to purchase in support of a 

new system depends, of course, not only on distribution of demand but also on the 

expected increase in cost for those spares if purchased at a later date. In an ideal 

situation having just the right number of spares over the lifetime of the program 

will bear the minimum costs. However, if costs are predicted to increase non- 

linearly over time then acquiring the optimal supply of spares becomes even more 

crucial. 

It is recognized that modelling non-linear costs as quadratic functions is only 
one approach. One might, for example, use cubic functions. A possible 
continuation of the work reported here is to explore optimal decisions for other 
non-linear cost relations. Another recommendation is to find optimal expressions 
for Q with results which would be more tractable than those found here. Another 
study continuing this work could address methods for solving for the coefficients 
C s | , C s 2 > etc.. 

This study has sought the optimal quantity of supply for both conditions of 
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risk and uncertainty about demand, when the cost function is quadratic. This 
situation has not been studied extensively in the past. Hopefully the results of this 
thesis will provide useful insights into the decision process while making 
affordability assessments of new acquisitions. 
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APPENDIX A 



To be determined is the maximum value of a maximization argument defined 
by an inequality. From page 39 we are given 

CJCn^CJQ') < C sl (,Q-+\?+CJQ-+l) 

and 

CJD ^ -<?*-!) < CJP m -Qf ♦ CJD^-Q-) . 

We wish to show that the maximum of the right hand side of 

MAX [ CJQ-f*CJ,Q') , C.fP m -Q'?*CJP„ a -Q’) } < 

MAX [ C SI (tQ‘+l?+C s2 (.Q-+l) , C 01 (,D mMX -Q*-l) 2 +C o2 (D aM1 -Q‘-l) ] , 

is 

C„(.Q-+lf+CJ.Q-+l) . 

Proof by contradiction: suppose 

C„(Q’+l?+CJQ-+l) < C.,(.D^-Q--lf.CJD^-<y-l) , 

then 
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C ol (.D^-Q'-Y?+CJP m -Q'-V) 

is greater than both 

CJQ-) 2 *CJQ-) and CJD^-Q'f+CJD^Q'), 
but this contradicts the one of the givens that 

is less than 

CJD m „-Q') 2 *CJD^-Q-) . 
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To be determined is the minimum value of a maximization argument defined 



by an inequality. From page 40 we are given 

CJQ'-Of + CJQ'-O) > CJQ'-lf+CJQ--!), 

and 

We wish to show that the minimum of the right hand side of 

UAX [ CJ.Q-f*CJ.Q-) , + } < 

is 

Proof by contradiction: suppose 

CJP^-Q'+lf + CJP^-Q'**) < CJQ‘-lY+CJQ--l) , 

then 
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is greater than both 



CJQ-) 2 *C s 2 (Q-) and C^-Of , 

but this contradicts one of the givens that 

C sjW“-l) 2 + C s2«?'-l) 



is less than 



CJQ‘f+CJQ-) 
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